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Principal component analysis of the reionization history from Planck 2015 data 

Wei-Ming Dai0 Zong-Kuan Guolll and Rong-Gen Ga0 
State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 

Chinese Academy of Sciences, P.O. Box 2735, Beijing 100190, China 

The simple assumption of an instantaneous reionization of the Universe may bias estimates of 
cosmological parameters. In this paper a model-independent principal component method for the 
reionization history is applied to give constraints on the cosmological parameters from recent Planck 
2015 data. We find that the Universe are not completely reionized at redshifts a > 8.5 at 95% 

CL. Both the reionization optical depth and the matter fluctuation amplitude are higher than 
but consistent with those obtained in the standard instantaneous reionization scheme. The high 
estimated value of the matter fluctuation amplitude strengthens the tension between Planck CMB 
observations and some astrophysical data, such as cluster counts and weak lensing. The tension can 
significantly be relieved if the neutrino masses are allowed to vary. Thanks to a high scalar spectral 
index, the low-scale spontaneously broken SUSY inflationary model can fit the data well, which is 
marginally disfavored at 95% CL in the Planck analysis. 

PACS numbers: 98.80.Cq 


1. INTRODUCTION 

Accurate measurements of the reionization history of 
the Universe plays an important role in understanding 
the early stages of structure and star formation. Fortu¬ 
nately, the evolution of the inter-galactic Lyman-alpha 
opacity measured in the spectra of quasars can provide 
valuable information on the reionization history [l| . The 
recent measurements suggest that the reionization of the 
inter-galactic medium was largely complete at redshift 
0 ~ 6 0- The rapid decline in the space density of 
Lyman-alpha emitting galaxies in the region z = 6—8 also 
implies a low-redshift reionization process Q. However, 
the detailed evolution of the electron ionization fraction 
with redshift is currently unknown due to the difficul¬ 
ties in probing the reionization state of the Universe at 
higher redshifts through the Lyman-alpha opacity. The 
most promising observation in the near future is detec¬ 
tion of the 21 cm transition of neutral hydrogen. The 21 
cm emission is a direct probe of the large-scale distribu¬ 
tion of neutral hydrogen during the epoch of reionization, 
and is, therefore, a complementary probe to the Lyman- 
alpha emission at lower redshifts. According to current 
understanding of reionization sources, high energy pho¬ 
tons may arise from the first star-forming galaxies and 
quasars. The uniqueness of these objects as the sources 
of photons to reionize the Universe is more challenging 
to establish. 

High precision measurements of the cosmic microwave 
background (CMB) anisotropies provide new insights 
into the reionization history. There are two main ef¬ 
fects of reionization on the CMB angular power spectra. 
The first effect produces a suppression of the acoustic 
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peaks in the CMB angular power spectra through the 
optical depth integrated over the whole reionization his¬ 
tory. Earlier reionization leads to the larger suppres¬ 
sion of the acoustic peaks. The second effect changes the 
shape of the large-scale polarization angular power spec¬ 
tra. This gives a characteristic bump in the spectra on 
scales larger than the horizon size at reionization. The 
position of the bump is proportional to the square root 
of the redshift at which the reionization occurs, while 
the amplitude is proportional to the optical depth Si- 
Therefore, measurements of the large-scale polarization 
angular power spectra can be used to constrain the reion¬ 
ization history ^ Q . 

Current constraint on the reionization optical depth 
from the Planck temperature data in combination with 
the low-multipole polarization data is r = 0.078 ± 0.019 
in the base ACDM model Q, which is smaller than r = 
0.089 ± 0.014 estimated by the WMAP 9-year data 
but consistent with that from WMAP polarization mea¬ 
surements cleaned for dust emission using 353 GHz po¬ 
larization maps from the high frequency instrument [lOj . 
These results are based on the assumption of an instan¬ 
taneous reionization. If such a simple assumption is not 
true, the estimated values of cosmological parameters are 
biased. Since there is currently no convincing model of 
the reionization history, it is important and necessary to 
constrain it in a relatively model-independent way. 

Hu and Holder proposed a principal component anal¬ 
ysis (PCA) of the reionization history to quantity the 
information contained in the large-scale A-mode polar¬ 
ization 0. The principal components are the eigenfunc¬ 
tions of the Fisher matrix that describes the dependence 
of the polarization angular power spectra on the reioniza¬ 
tion history. In practice, no more than the first five prin¬ 
cipal components are needed to parameterize the reion¬ 
ization history 0. This approach has been applied to 
WMAP data and simulated future data [l^ - [l^ . It is 
recently claimed that the Harrison-Zel’dovich primordial 
spectrum is consistent with the WMAP 7-year data in 
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such a general reionization scenario [3, [l^ . In this pa¬ 
per, we perform the principal component analysis of the 
reionization history from recently released Planck 2015 
data and study the impact of the reionization history 
on the estimates of the cosmological parameters. We 
also investigate the bounds on the neutrino masses and 
primordial tensor perturbations in a general reionization 
scheme. 

The paper is organized as follows. In Section [2] we 
describe the principal component method for parameter¬ 
izing the reionization history that we will adapt in our 
analysis. The results are presented in Section [3l The last 
section is devoted to our conclusions. 


2. MODEL 



Since there is no convincing physical model of the 
reionization history, the instantaneous reionization his¬ 
tory is usually assumed when investigating cosmologi¬ 
cal constraints from CMB measurements. Such a one- 
parameter phenomenological model may bias the esti¬ 
mates of the optical depth and other cosmological pa¬ 
rameters. In general, the hydrogen ionization fraction 
Xe{z) is a function of redshift. Binning the hydrogen ion¬ 
ization fraction in redshift bins is a good idea to parame¬ 
terize the reionization history. However, the degeneracy 
between the binned ionization fraction leads to weak con¬ 
straints [ 2 ^. The principal component analysis provides 
a solution to this problem Q. 

The PCA approach is often adopted to reduce the num¬ 
ber of variables while containing most of information of 
the system. It converts a set of correlated variables into 
a set of linear uncorrelated variables by an orthogonal 
transformation. The principal components are the eigen¬ 
vectors of the diagonalized matrix, which are picked out 
according to the corresponding eigenvalues. Intuitively, 
the larger the eigenvalue is, the more information the 
corresponding eigenvector gives. 

Following ,we parameterize the reionization his¬ 
tory. Consider a binned ionization fraction Xe(zi), i S 
{1,2,..., iV^}, with redshift bins of width Az = 0.25 
spanning ^ z < Zmax- The bin width of Az = 0.25 
is sufficiently small so that the final results we obtain 
are independent of the redshift binning m- Here we 
take Zmin = 6 and Zmax = 30. With the definition 
Zl = Zmin -I- Az and Zn^ = .^max — Az so that A 2 -I- 1 = 
(^max — Zmin)/Az. For z > Zmaxj the ionization frac¬ 
tion is set to the residual fraction of recombination, while 
Xe = 1.0 for 3 < z < Zmin Consistent with astrophysical 
observations 0 and Xe = 1.08 for z < 3 when the Helium 
reionization is taken into account. The principal compo¬ 
nents of Xe{zi) are the eigenfunctions of the following 
Fisher matrix Fij, 


/ i\ dlnCrdlnCr 

^ V 27 dXe{Zi) dXe(Zj) 


( 1 ) 


FIG. 1: Fractional power spectrum response to a delta func¬ 
tion perturbation of unit amplitude at 6.25 < Zi < 29.75. 
Here we have chosen Az = 0.25. 


which describes the dependence of the polarization spec¬ 
trum on the ionization fraction Xe(zi). To calculate 
the Fisher matrix, we should set a fiducial reionization 
history Xg^(zi) as well as other cosmological parameters. 
The fiducial model parameters are chosen as the esti¬ 
mated values of cosmological parameters from the Planck 
2015 data including the temperature and polarization 
power spectra for the base ACDM model Q, listed in 
Table m We adopt Xg‘^(zi) = 0.1 in the range of redshifts 
6 < z < 30 to get the same optical depth as the instanta- 
nous reionization. Actually our results are not sensitive 
to the fiducial model. Since the effect of Xe{zi) on 
of can be neglected on high i, as shown in Fig. (TJ we 
can cut off £ safely to imax = 100 in the Fisher matrix. 
The perturbation to the fiducial x^ is a delta function 
correlated with the nearby region 0. 


TABLE I: Fiducial model consistent with results from the 
2015 Planck temperature and polarization spectra for the 
standard ACDM model. 




h r 

AsC ris 

0.02225 

0.1198 

0.6727 0.079 

1.882 X lO'^* 0.9645 


The Fisher matrix Fij is decomposed as 


Fij = {Nz -|- 1 ) S^{zi)<j^ S^{zj ), ( 2 ) 

11=1 

where are the inverse eigenvalues and S^{z) are the 
eigenfunctions that satisfy the orthogonality and com¬ 
pleteness relations 


L 


min 

N, 


dzSi_i(^z)Si/(^z') — (-S^max ; (3) 

( 4 ) 


Sfi{zi)Sfx{zj) — {Nzl)Sij. 
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So these eigenfunctions form a set of orthogonal and com¬ 
plete bases. An arbitrary reionization history can be de¬ 
composed into a sum of the eigenmodes, 

a;e(z) =x“(z) + ^m^S'^(z), (5) 

where are the amplitudes of the principal compo¬ 
nents for a particular reionization history. We order the 
eigenmodes so that the smallest has the lowest index 
number starting at /i = 1. Since the effect of each 
eigenmode on becomes smaller as fi increases, we 
shall pick out the first few eigenmodes. It is enough to 
choose the first 3 — 5 eigenmodes shown in Fig. [2] because 
of the existence of the cosmic variance 

ACf® 

^EE 

so that the error in Cf^ from high /i eigenmodes is 
smaller than the cosmic variace. 


2£-h 1 ’ 


( 6 ) 


S,(z) 



FIG. 2: The five lowest-variance principal components of 
Xe{z) over the redshift range of 6.25 < z < 29.75. The vari¬ 
ances increase with the index jj,. 


In our analysis, we parameterize the reionizaton his¬ 
tory as with the five amplitudes for /r = 1,..., 5. 
We consider a spatially-flat ACDM model described by 
a set of cosmological parameters 

0s, As, ris}, 

where and are the physical baryon and cold 

dark matter densities relative to the critical density, 0 s 
is the ratio of the sound horizon to the angular diameter 
distance at decoupling, Ag is the amplitude of the pri¬ 
mordial power spectrum of scalar perturbations and ris is 
the scalar spectral index. Moreover, we consider two one- 
parameter extensions to the ACDM model, ACDM-|-Em^ 
and ACDM-|-r, where is the total mass of neutri¬ 
nos and r is the tensor-to-scalar ratio. When the tensor 
contribution to the CMB angular power spectra is con¬ 
sidered, the inflationary consistency relation, rit = —r/8, 
is assumed where rit is the tensor spectral index. For 
comparison with the Planck results, r is defined at the 
pivot scale fc» = 0.002 Mpc“^ (denoted ro,oo 2 )- 


We modified the Boltzmann CAMB code [2l| to ap¬ 
propriately incorporate the PCA reionization history and 
used the publicly available CosmoMC package to explore 
the parameter space by means of Monte Carlo Markov 
chains technique 

In our analysis we use recently released Planck like¬ 
lihood code and data, including the Planck low-Z like¬ 
lihood at multipoles 2 < I < 29 and high-Z Plik likeli¬ 
hood at multipoles I > 30 based on pseudo-C; estima¬ 
tors. The former uses foreground-cleaned LFI 70 GHz 
polarization maps together with the temperature map 
obtained from the Planck 30 to 353 GHz channels by the 
Gommander component separation algorithm over 94% 
of the sky. The polarization part of this likelihood is 
denoted as “lowP”. The latter uses 100 GHz, 143 GHz, 
and 217 GHz half-mission cross-power spectra, avoiding 
the Galactic plane as well as the brightest point sources 
and the regions where the CO emission is the strongest. 
“Planck TT-|-lowP” denotes the combination of the TT 
likelihood at Z > 30 and a low-Z temperature-polarization 
likelihood. “Planck TT,TE,EE-|-lowP” denotes the com¬ 
bination of the likelihood at Z > 30 using TT, TE, and 
EE spectra and the low-Z temperature-polarization like¬ 
lihood. 


3. RESULTS 

Table m summarizes the constraints on the reioniza¬ 
tion history and cosmological parameters from the Planck 
TT-|-lowP and Planck TT,TE,EE-|-lowP data. The re¬ 
sults from the temperature and polarization power spec¬ 
tra are consistent with those from Planck TT-|-lowP, 
but with increased precision. Since temperature-to- 
polarization leakage corrections are not considered in the 
2015 Planck TE and EE power spectra Q, in our analysis 
Planck TT-t-lowP is usually quoted. Actually the reion¬ 
ization history is not sensitive to the TE and EE power 
spectra at Z > 30. 

The mean posterior values and 68 % confidence limits 
of are listed in Table im Compared to WMAP data, 
Planck data place strong constraints on the amplitudes 
of the principal components. The marginalized 2D con¬ 
tours and posterior probability distributions are plotted 
in Eigure[71 As expected, the correlations between {m^} 
are weak and the probability distributions are Gaussian. 
Moreover, we reconstruct the reionization history from 
Monte Carlo chains using the PCA approach in Figure |31 
The blue curve is the mean value of the hydrogen ioniza¬ 
tion fraction Xe{z) in the redshift range 6.25 < z < 29.75, 
with 68 % (dark color) and 95% (light color) confidence 
regions. As we can see the Universe is not completely 
reionized at z > 8.5 at 95% CL. This result is consis¬ 
tent with the estimated value of the reionization redshift 
Zre = 8 . 8 ti ;4 derived from Planck TT-|-lowP-|-lensing in 
the instantaneous reionization model. 

The constraint on the optical depth of the PCA reion¬ 
ization is r = 0.092 ±0.017, lager than that of the instan- 
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TABLE II: 68% constraints on the PCA reionization and other cosmological parameters from 2015 Planck CMB angular power 
spectra. For the total mass of neutrinos and tensor-to-scalar ratio, the 95% upper limits are given. 


Parameters 

ACDM 

ACDM 

ACDM+r 

ACDM + Smj/ 

Planck TT+lowP 

Planck TT,TE,EE+lowP 

Planck TT + lowP 

Planck TT + lowP 


0.02223 ± 0.00023 

0.02224 ± 0.00016 

0.02225 ± 0.00023 

0.02210 ± 0.00030 


0.1192±0.0021 

0.1196 ± 0.0015 

0.1190 ± 0.0021 

0.1201 ± 0.0024 

1006»mc 

1.04092±0.00047 

1.04079 d= 0.00032 

1.04095 ± 0.00047 

1.04066 ± 0.00056 

T 

0.092±0.017 

0.092 d= 0.015 

0.092 ± 0.017 

0.095 ± 0.017 

Tig 

0.9689±0.0064 

0.9663 d= 0.0049 

0.9698 ± 0.0064 

0.9662 ± 0.0075 

< 

o 

o 

3.115±0.033 

3.117 ± 0.030 

3.115 ± 0.032 

3.121 ± 0.034 

Hq (km a ^ Mpe ^) 

67.54±0.96 

67.34 ± 0.65 

67.64 ± 0.96 

64.85 ih 3.34 

^8 

0.838 ± 0.014 

0.840 ± 0.012 

0.839 ± 0.014 

0.788 ± 0.063 


0.3116 ± 0.0132 

0.3144 ± 0.0091 

0.3104 ± 0.0131 

0.3507 ± 0.0525 

’'0.002 

- 

- 

< 0.091 

- 

Em,, (eV) 

- 

- 

- 

< 1.03 

mi 

U.U5by_Q Q534 

U.U485_q Q5^Q 

0.0528j:“;“®“| 

n •0712 

0.0663_q Q5g3 

m2 

-0 0457+°'^^^^ 
‘^•^^°'_0.1091 

U.Ub8y_Q 

-0.05431®;“?! 

n .1300 

-0.0536_q ^ggg 

m3 

u.u^yo_Q 

U.U4uy_Q 

0 0586”^*^■ 
U.U5Sb_Q 

U.Ub42_Q j^g^3 

m4 

-0 010+°'^®°3 
U.UiU_Q 

U.UiyU_Q ^ggg 

U.Ui2<-_Q ^^27 


ms 

'^■‘^^^^-0.1570 

U.U4Ub_Q ^44g 

U.U519_q j^g^g 

n nccn”f"0 • 1688 

0.0560_q j^ggg 


taneous reionization obtained in the 2015 Planck analy¬ 
sis, but with nearly the same uncertainty Q. The reason 
is that the ionization fraction at high redshift makes a 
larger contribution to the optical depth than at low red- 
shift. The optical depth between any two redshifts zi 
and Z 2 is an integration of Xe, 


t{zi,Z2) = 0.0691(1 



dz 


(1 + ^)^ 

H{z)/Ho 


Xe{z) 


(7) 

where Yp is the helium abundance. Since the Hubble pa¬ 
rameter H falls off as (l + z)^^'^ in the matter domination 
epoch, the same Xe can contribute more to the optical 
depth at higher redshift. It is known that the amplitude 
of primordial spectrum of scalar perturbations degen¬ 
erates with optical depth r in the form Ase~‘^'^ on small 
scale [l^, which means that a large t leads to a large 
Ag. Hence a higher value of the matter fluctuation am¬ 
plitude, (Tg = 0.838 ±0.014, is favored in the PCA reion¬ 
ization than in the instantaneous reionization, as shown 
in the left panel of Figure ID This means that such a 
general reionization history strengthens the tension be¬ 
tween Planck CMB observations and some astrophysical 
data, such as cluster counts and weak lensing. As shown 
in the right panel of Figure IH this tension can signifi¬ 
cantly be relieved if the neutrino masses are allowed to 
vary, thanks to the degeneracy between the reionization 
history and the neutrino masses. Larger neutrino masses 
lead to a lower erg through the effects of neutrino free 
streaming on structure formation. The constraint on the 
total mass of neutrinos can be relaxed to Yrrip < 1.03 
eV at 95% CL in the PCA reionization scenario. In Fig¬ 
ure [5] we also show the correlation of the neutrino masses 
with the Hubble constant. The constraint has a broad 
tail to high masses, which illustrates the acoustic scale 
degeneracy with Hq. 



Z 


FIG. 3: Reconstructed reionization history with the PCA 
method from PlanckTT-flowP. The blue curve is the mean 
value of the hydrogen ionization fraction Xe[z) in the redshift 
range 6.25 < z < 29.75, with 68% (dark color) and 95% (light 
color) confidence regions. 


It is claimed in [T^ that the Harrison-Zel’dovich 
primordial spectrum (i.e., scale-invariant spectrum) of 
scalar perturbations is consistent with the WMAP 7-year 
data in the PCA reionization scenario. However, our 
analysis indicates that the Harrison-Zel’dovich spectrum 
is disfavored by Planck TT±lowP at 4.9 a CL, although 
the scalar spectral index is shifted towards slightly higher 
values. To test inflationary models the contribution 
of primordial tensor fluctuations to the CMB angular 
power spectra is taken into account. Figure |6] depicts the 
marginalized 2D contour in the ro.oo 2 — ng plane, with 
68% and 95% CL, for the PCA (blue) and instantaneous 
(red) reionization history, derived from Planck TT±lowP. 
We find that the allowed contour is shifted towards to 
higher values of rig, compared to the instantaneous reion- 
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FIG. 4: Marginalized 2D contours in the as — fim plane, with 68% and 95% CL, for the ACDM model (left panel) and 
ACDM-|-Emiy (right panel) with the PCA (blue) and instantaneous (red) reionization history, derived from Planck TT+lowP. 




FIG. 5: Marginalized 2D contour in the ILo —plane, with 
68% and 95% CL, for the ACDM-I-Em,/ model with the PCA 
(blue) and instantaneous (red) reionization history, derived 
from Planck TT+lowP. 


ization. This shift is important for the constraints on in¬ 
flationary models. For example, it is pointed out in 
that the class of inflationary models with a power-law po¬ 
tential for p <C 1 and the low-scale spontaneously bro¬ 
ken SUSY inflationary model are marginally disfavored 
by Planck data at 95% CL. In the PCA reionization sce¬ 
nario, however, these models are in agreement with re¬ 
cent CMB observations. Moreover, the constraints on 
the tensor-to-scalar ratio is tighter compared to the in¬ 
stantaneous reionization because the large amplitude of 
scalar spectrum suppresses the tensor-to-scalar ratio. 


FIG. 6: Marginalized 2D contour in the ro.oo 2 — n-s plane, with 
68% and 95% CL, for the ACDM-|-r model with the PCA 
(blue) and instantaneous (red) reionization history, derived 
from Planck TT-|-lowP. 


4. CONCLUSIONS 

We have considered the principal component 
parametrization of the reionization history and studied 
the effects of the PCA reionization on the estimates of 
cosmological parameters using recently released Planck 
data. From the reconstructed reionization history we 
found that the Universe is not completely reionized 
at z > 8.5 at 95% CL. This conclusion can be tested 
by future astrophysical experiments such as SKA. The 
reionization optical depth is higher than but consistent 
with that obtained in the instantaneous reionization 
scheme. The higher value of r leads to the higher 
value of As, hence the higher erg, which strengthens the 
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FIG. 7: Marginalized 2D contours (68% and 95% CL) and posterior distributions for the amplitudes of the principal components 
derived from Planck TT+lowP. 


tension between Planck CMB observations and some 
astrophysical data, such as cluster counts and weak 
lensing. But the tension can significantly be relieved if 
the neutrino masses are allowed to vary. The constraint 
on the total mass of neutrinos from CMB data is relaxed 
due to the degeneracy with the reionization history. 
Thanks to the shift of the scalar spectral index to higher 
values in the PCA approach, the low-scale spontaneously 
broken SUSY inflationary model, which is marginally 
disfavored at 95% confidence level in the Planck analysis, 
can fit the data well. 
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